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ABSTRACT 

Smaller terrestrial planets (< 0.3 M®) are less likely to retain the substantial atmospheres and 
ongoing tectonic activity probably required to support life. A key element in determining if sufficiently 
massive "sustainably habitable" planets can form is the availability of solid planet-forming material. 
We use dynamical simulations of terrestrial planet formation from planetary embryos and simple 
scaling arguments to explore the implications of correlations between terrestrial planet mass, disk 
mass, and the mass of the parent star. We assume that the protoplanetary disk mass scales with 
stellar mass as M^sk ^ f where / measures the relative disk mass, and 1/2 < h < 2, so that disk 
mass decreases with decreasing stellar mass. We consider systems without Jovian planets, based on 
current models and observations for M stars. We assume the mass of a planet formed in some annulus 
of a disk with given parameters is proportional to the disk mass in that annulus, and show with a 
suite of simulations of late-stage accretion that the adopted prescription is surprisingly accurate. Our 
results suggest that the fraction of systems with sufficient disk mass to form > 0.3M ffi habitable 
planets decreases for low-mass stars for every realistic combination of parameters. This "habitable 
fraction" is small for stellar masses below a mass in the interval 0.5 to 0.8 M©, depending on disk 
parameters, an interval that excludes most M stars. Radial mixing and therefore water delivery are 
inefficient in lower-mass disks commonly found around low-mass stars, such that terrestrial planets 
in the habitable zones of most low-mass stars are likely to be small and dry. 

Subject headings: astrobiology — planetary systems: formation — planetary systems: protoplanetary 
disks — stars: low-mass, brown dwarfs 



1. INTRODUCTION 

The circumstellar habitable zone (HZ) is the distance 
annulus within which a terrestrial planet can maintain 
liquid water on its surface. The average distance of the 
HZ from the parent star depends primarily on the star's 
luminosity, and is closer-in and narrower around lower- 
mass stars (Kasting et al. 1993). For a planet to remain 
habitable, it must control its surface temperature over 
long timescales, possibly via the C02-carbonate cycle 
(Walker et al. 1981), which is enabled by a relatively 
thick atmosphere, a hydrological cycle, and active plate 
tectonics. These planetary characteristics are less likely 
to be sustained over long periods by less massive terres- 
trial planets, or those with low water contents. Whether 
potentially habitable planets of sufficient mass will form 
depends on the surface density of rocky material in their 
parent protoplanetary disks. Disks with higher surface 
densities tend to form a smaller number of more massive 
planets (Wetherill 1996; Raymond et al. 2004, Kokubo 
et aL2006), and can increase terrestrial planet water 
content by scattering more distant, water-rich material 
into the HZ (Morbidelli et al. 2000; Raymond et al. 2004, 
2005a). 

Here we investigate the formation, masses and water 
contents of terrestrial planets as a function of stellar mass 
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by assuming a parameterized relationship between disk 
and stellar mass. We adopt a simple model relating the 
HZ planet mass to the disk mass, based on results in 
the literature. We use a set of numerical simulations of 
terrestrial planet growth around stars between 0.2 and 1 
M© to test the validity of the model. The simulations 
agree well with the model, and justify an exploration 
of the consequences of such a planet-disk-stellar mass 
relation, and in addition allow us to consider how water 
delivery is affected as a function of stellar mass. 

Our main result is to derive probabilistic limits on the 
stellar mass below which sustainably habitable planets 
cannot form because of insufficient disk mass available 
for accretion to a critical planet mass of 0.3 M®. In ad- 
dition, we show how this stellar mass limit depends on 
the parameters of the disk-star mass relationship, which 
is currently poorly constrained. In a statistical frame- 
work, we show that the fraction of stars capable of form- 
ing habitable planets in situ decreases quickly with de- 
creasing stellar mass, and that terrestrial planets in the 
habitable zones of M stars (M* < 0.5 M©; e.g., Reid & 
Hawley 2000) are likely to be small and dry. An interest- 
ing test of our models comes from the recent detection 
of a ~ 5 M© planet close to the habitable zone of the 
M dwarf Gliese 581 (Udry et al. 2007). The implications 
of this system for our formation models are discussed in 
section 4. 

To simplify our calculations, we assume there are no 
gas giant (Jovian) planets present, and we test the effects 
of this assumption in section 3.3. Systems without giant 
planets may represent the majority of planetary systems 
with habitable planets. Available estimates of the frac- 
tion of solar-like stars with massive planets (Tabachnick 
&; Tremaine 2002, Lineweaver & Grether 2004, Beer et 
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al. 2004; Fischer & Valenti 2005) allow for 70-90% of stars 
to have planetary systems without Jovians, and an even 
larger fraction for lower-mass stars (Endl et al. 2006; But- 
ler et al. 2006). In addition, Scholz et al. (2006) point out 
that estimated submm and mm outer disk masses show 
that only a small fraction of low-mass stars or substel- 
lar objects have sufficient disk mass to form a planet 
with the mass of Jupiter (but see Hartmann et al. 2006). 
The lack of a correlation between the presence of debris 
disks and the metallicity of the parent star (Greaves et 
al. 2006; Beichman et al. 2006; Moro-Martm et al. 2007), 
unlike the correlation known for exogiant planets (Gon- 
zalez 1997, Santos et al. , Fischer & Valenti 2005), also 
supports the likelihood of gas-giant-free planetary sys- 
tems. There are observational (Endl et al. 2006; Butler 
et al. 2006) and theoretical (Ward & Hahn 1995; Adams 
et al. 2004; Laughlin et al. 2004, Ida & Lin 2005) reasons 
to think that the fraction of M stars with Jovian planets 
is significantly smaller than for solar-mass stars. 

2. METHODS 

2.1. Habitable Planet Mass Limit 

Following Williams et al. (1997), we assume a lower 
mass limit for planetary habitability that supports ac- 
tive plate tectonics for several Gyr. Significant mass is 
a necessary, but not sufficient, requirement for sustain- 
able planetary habitability (e.g., Lissauer 1999). Using 
the radioactive flux limit from Williams et al. (1997), the 
critical planetary mass scales like p~ 2 e 3At , where p is the 
planet's bulk density, A — 1.5 x I0 10 yr~ 1 for 238 U, and t 
is the duration of tectonic activity. For an active lifetime 
of 5 Gyr, and iterating on the density to be consistent 
with the resulting critical mass, we get a critical mass of 
0.3 M ffi for a density of 4.5 gem' 3 . This result is sen- 
sitive to the plate tectonic activity timescale associated 
with sustained habitability, which is fairly arbitrary. The 
average age of stars in the Galactic disk is about 5 Gyr, 
so the choice of 5 Gyr allows about half of such stars 
with sustainably habitable planets to be on the main se- 
quence today, for a roughly constant star formation rate. 
Throughout the paper we consider 0.3 M ffl as our critical 
mass limit for planetary habitability, and we discuss the 
implications of changing this limit in section 4. 

2.2. Protoplanetary Disk Properties 

We assume the surface density of protoplanetary disks, 
S, to scale with heliocentric distance r as: 




where £i is the surface density of solids at 1 AU in 
the minimum-mass solar nebula (MMSN) model (ss 
6g cm' 2 ), ft adjusts the density scaling with M+, / is a 
scaling factor for the disk (/=1 for MMSN), Z is the stel- 
lar metallicity (1.0=solar), and a moderates the steep- 
ness of the density profile (a = 3/2 in the MMSN model; 
Hayashi 1981, but more recent reconstructions using an 
accretion disk model as a constraint gives a — 1/2; Davis 
2005). Most disk models that include viscous and irra- 
diation heating give a between about 1/2 and 1 for the 
inner disk (e.g. D'Alcssio et al. 1998; Garaud & Lin 2007, 
Dullcmond et al. 2007), and high- resolution submillimc- 
ter observations suggest a between 1/2 and 1 in the outer 



disk (Mundy et al. 2000, Looney et al. 2003, Andrews 
and Williams 2006). We take a = 1 as a fiducial value 
but consider values in the range 1/2 to 3/2. 

Present estimates of disk masses in nearby star-forming 
regions (see Andre & Montmerle 1994, Eisner & Carpen- 
ter 2003, Andrews & Williams 2005) give a large spread, 
with median or average masses usually somewhat below 
the MMSN value, but a significant fraction more massive, 
including a very massive disks observed in both Orion 
and Taurus (e.g. Williams et al. 2005, Eisner & Carpen- 
ter 2006). At the other end of the mass spectrum, the 
estimates by Williams et al. (2005) of mean masses for 18 
Orion proplyds below the 3 — a detection limit yielded 
8 x 10~ 4 M Q , or only 3 M e in solids for a gas-to-dust ratio 
of 100. There also exists a factor of 3-5 range in metal- 
licity for thin-disk stars (Nordstrom et al. 2004), which 
adjusts the effective gas-to-dust ratio. 

There is considerable uncertainty in the scaling be- 
tween stellar mass M* and disk mass M<n s k- The most 
promising lead comes from estimates of outer disk masses 
using submm and mm techniques (Beckwith et al. 1990), 
which remain uncertain because of adopted dust tem- 
peratures, opacities, distances, and radiative transfer 
model. Scholz et al. (2006) estimated disk masses for 
brown dwarf disks and combined their results with previ- 
ous determinations for stellar and substellar objects from 
0.02 to 3 Mq. Their Fig. 3 shows that the ratio of disk 
mass to stellar mass is roughly constant with stellar mass, 
with a large amount of scatter. This hints that ft ~ 1, 
with a large range in /. Many of the data points from 
Scholz et al. (2006) are only upper limits, but we think 
extreme cases ft = or ft > 2 would be apparent in this 
data. 

Additional clues come from the observed stellar ac- 
cretion rate M, which scales as M oc M 2 , with a 
large dispersion (Muzcrollc et al. 2003, 2005, Natta et 
al. 2006). It might therefore seem natural to assume 
that Mdisk oc M 2 , i.e. ft=2, assuming that the viscous 
timescale is independent of stellar mass (e.g., Ida & Lin 

2005) . On the contrary, nearly all accretion disk models 

predict M oc with a <~ 1. This matter is still unre- 
solved; several solutions have been proposed (see Padoan 
et al. 2006; Alexander and Armitage 2006; Dullcmond et 
al. 2006; Hartmann et al. 2006; Gregory et al. 2007) in- 
cluding selection/detection limitations (Clarke & Pringle 

2006) . Most of these models are consistent with ft « 1. 
Although ft is currently poorly constrained, we empha- 
size its importance for understanding whether sustain- 
able HZ planets can form around low-mass stars. Given 
the results of Scholz et al. (2006), we choose ft = 1 as our 
fiducial case, but explore the effects of ft from 0.5 to 2. 

2.3. Simulations and Models 

We performed dynamical simulations to evaluate the 
scaling between terrestrial planet mass and disk prop- 
erties. For all runs, our starting conditions reflect our 
fiducial case, with fZ = 1, ft = 1, and a = 1. We 
vary the mass of the central star from 0.2 to 1 M . 
Given the surface density profile from Eq (1), we assume 
that planetary embryos have formed throughout the in- 
ner disk of each star. We assume that the disk evolved 
following the standard model of terrestrial planet growth 
(e.g. Chambers 2004): grains coalesced to form km- 
sized planetesimals; and embryos formed from planetes- 
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imals via runaway and oligarchic growth (e.g., Kokubo 
& Ida 1998), spaced by 5-10 mutual Hill radii (Rh,™ = 
0.5(ai + a 2 )(M 1 +M 2 /3M lt ) 1 / 3 ; oi and d2 are the orbital 
radii and Mi and Mi the masses of two adjacent em- 
bryos). For the case of 1 M , we generated a population 
of ~ 75 embryos from 0.5 to 4 AU, totaling 4.95 M®. 6 
For other stellar masses, we assume the population of 
embryos has the same temperature, as we are interested 
in HZ planets (defined to have the same temperature), 
and the feeding zone boundaries should scale in the same 
way as the HZ mean distance. We therefore scale the in- 
ner and outer boundaries by the stellar flux (i.e., as the 

1/2 

luminosity Lj ), using a mass- luminosity relation from 
Scalo et al. (2007) that is a fit to data from Hillenbrand 
& White (2004): y = 4.101a; 3 + 8.162x 2 + 7.108a; + 0.065, 
where y = loglL*) and x = log(M- k ). Table 1 summa- 
rizes our starting conditions for each stellar mass. Note 
the variations in the number of embryos N included in 
each simulation, from ~ 75 for the 1 M simulations to 
almost 200 for the 0.2 M simulations. This variation 
in N is due to the Hill radius being very small at the 
small orbital distances studied for lower-mass stars. In 
other words, models of embryo growth predict that there 
really are more embryos in these inner disks than in the 
more distant regions studied for higher-mass stars (e.g., 
Kokubo & Ida 1998). Kokubo et al. (2006) showed that 
the bulk properties of accreted planets are not particu- 
larly sensitive to such variations in TV, so we are not con- 
cerned that our choice of embryo formation models will 
affect our results. Note also that the 0.2 M star has less 
than 0.3 M in its entire disk, so it serves primarily to 
tell us about the efficiency with which the disk mass is 
used to make planets under these conditions, but cannot 
form a habitable planet according to our criterion. Each 
numerical experiment was run independently four times 
starting from different random initial embryo positions 
for a total of twenty simulations. We give embryos ran- 
dom starting eccentricities <0.02 and inclinations <0.1 
degrees. 

Embryos are assigned water contents based on values 
from our Solar System, where asteroids beyond ^2.5 AU 
contain significant quantities of water (Abe et al. 2000; 
see Fig. 2 in Raymond et al. 2004). We assumed that this 
boundary between "dry" and "wet" embryos (called the 
"water line" in Table 1) scales with the stellar flux. 7 Em- 
bryos beyond this boundary contain 5% water by mass, 
and those inside are dry (Table 1). We integrated each 
simulation for 200 Myr using the hybrid integrator Mer- 
cury (Chambers 1999). Timesteps were chosen to sample 
the innermost body's orbit at least twenty times (e.g., 
Rauch & Holman 1999; Levison & Duncan 2000; see Ta- 
ble 1), and so varied with each set of simulations, from 

6 Note that the outer edge of 4 AU for the embryo disk is chosen 
based on previous simulations showing that the feeding zone of a 
terrestrial planet in the habitable zone does not extend beyond 4 
AU (Raymond et al. 2007). 

7 We use the term "water line" instead of the more common 
"snow line" because we are not necessarily assuming this to be the 
location where the temperature drops below ~ 170 K and water ice 
can condense (2.7 AU in the MMSN; Hayashi 1981). Rather, we 
are simply assuming this boundary to divide dry and wet material. 
The water line may actually be located somewhat interior to the 
snow line, as water-rich bodies can be shifted inward via either gas 
drag or by eccentricity pumping during embryo formation (e.g., 
Cyr et al. 1998; Kokubo & Ida 1998). 



6 days for 1 M Q to 0.2 days for 0.2 M . Collisions were 
treated as inelastic mergers conserving water. 

3. RESULTS 
3.1. Terrestrial Planet Mass vs Stellar Mass 

Figure [T] shows the evolution of a simulation for a 
0.6 M star: the disk is excited by gravitational pertur- 
bations among the embryos. As eccentricities increase, 
orbits cross and collisions occur. In time, a few planets 
grow and the number of bodies dwindles. Embryos may 
be scattered far from their original locations, sometimes 
delivering water-rich material to planets in the inner re- 
gions. Water delivery occurs relatively late, because mul- 
tiple scattering events are needed for significant radial 
movement (Raymond et al. 2007). In this case some wa- 
ter delivery occurred: a 0.21 M planet formed at 0.41 
AU, just beyond the outer boundary of the HZ, accreted 
two water-rich embryos originating beyond the water line 
at 0.61 AU. However, no water was delivered to any plan- 
ets in the HZ in this case. This is a typical outcome for 
the low-mass disks expected to be common around low- 
mass stars (see below). 

Figure O shows the final outcome of ten simulations, 
with the Solar System included for scale (the Earth's 
water content is ~ 10~ 3 bymass; Lecuyer et al. 1998). It 
is clear that, for our assumptions of h = 1, a = 1 and 
h = 1, terrestrial planets are much smaller around low- 
mass stars. In addition, planets that form in the HZs of 
low-mass stars tend to be dry, and more closely spaced. 
Note that, although we do not follow their orbits for a 
full 5 Gyr, we assume water-rich planets that form in the 
HZ to be potentially habitable. Very late-stage instabili- 
ties were not seen in any of these simulations, nor in the 
>Gyr integrations from Raymond et al. (2005a, 2006a), 
but the potential disruption of the system at times as late 
as 5 Gyr cannot be ruled out. In addition, a late-stage 
instability could potentially alter the orbit of a distant 
planet, causing it to collide with a planet in the HZ. Such 
an event could even deliver water at a very late stage. Al- 
though this is certainly possible, we have found this type 
of event to be rare in previous long-term simulations. In 
addition, the source of such a rogue planet would have to 
be quite far out, because accretion tends to occur faster 
closer to the star, and crossing orbits are needed to cause 
a strong scattering event. 

Figure [3] shows the mean mass of simulated planets 
that formed in the HZ as a function of stellar mass, with 
error bars representing the range of outcomes from the 
four simulations for a given stellar mass M*. The solid 
curve shows the prediction of a simple scenario in which 
the mass of a planet in the HZ is proportional to the 
mass contained within the HZ annulus, such that 

M , 3C El ? ZM * (V-a _ r 2-c \ (2) 
wipl IX 2- a \ HZ > out HZ,in) ' ^ ' 

where rnz,in and rRZ,out are the inner and outer bound- 
aries of the HZ (see Table 1). This model is a very 
simple, but not unreasonable, approximation for planet 
mass, and the quality of the fit is remarkable - note that 
this model was calibrated such that a disk with JZ = 1 
will form a 1 M® planet in the HZ around a solar-mass 
star. Kokubo et al. (2006) showed that the planet mass 
scales roughly linearly with the available mass, albeit 
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for a fixed stellar mass, and we independently found 
the same result for M* = 0.4 M©. In fact, Kokubo et 
al. (2006) found a slightly stronger than linear correla- 
tion, M p i oc M°fJr 1A . Fig [2] shows that the mean in- 
terplanetary spacing decreases somewhat for lower-mass 
stars, and that there are slightly more planets for the 
lower-mass stars. We know that the total mass in the 
habitable zone, Mjjz, is equal to the number of planets, 
N, times the average planet mass, M p i. Our simulations 
suggest that N oc or so. Indeed, a model with 

M p i oc M}jz provides a fit that is comparable to the one 
in Fig. [21 However, for the remainder of the paper we 
assume that M p i oc Mhz (Eqn 2). The reason for this 
assumption has to do with the goals of the paper. We 
are attempting to constrain the locations in — h — fZ 
parameter space that might harbor potentially habitable 
planets with M p i > 0.3 M®. To be conservative in our 
evaluations, we prefer to slightly overestimate, rather 
than underestimate, M v \. The difference between the 
two estimates is negligible for larger stellar masses, but 
is as much as ~ 40% below 0.1 M©. 

The planet mass decreases monotonically with stellar 
mass for all reasonable parameter values (Eq. 2; Fig [3]), 
with a scatter in the details of a given system based 
on the stochastic nature of the accretion process (e.g., 
Wetherill 1996). Only for very steep density profiles 
(a > 2) or reversed disk mass scalings (h < 0) can the 
planet mass increase at lower stellar masses. These ef- 
fects are the result of the strong dependence of the HZ's 
location's on stellar luminosity, and therefore on stellar 
mass. For a given value of /, Z, a, and h, there ex- 
ists a stellar mass limit below which the formation of a 

> 0.3 M© planet in the HZ is unlikely. For a = 1 and 
h = 1, this limit ranges from 1 M© for fZ < 0.3 to 0.74 
M© for fZ = 1 to 0.43 M© for fZ = 5. These limits 
clearly depend on the critical mass for habitability; for 
instance, the limit is 0.53 M© for the fZ = 1 case if the 
critical habitable mass is 0.1 M®. Recall that fZ repre- 
sents a scaling of the disk mass, i.e., the disk's relative 
mass / times the relative abundance of solids, assumed 
to scale with the stellar metallicity Z. 

Figure [4] shows the location in M* — fZ space where 
planets > 0.3 M© can form in the HZ, assuming a = 1. 
Each curve corresponds to a given value of h; planets 

> 0.3 M® form above and to the right of each curve. 
More massive or metal-rich disks can form habitable 
planets around lower-mass stars. In addition, it is easier 
to form > 0.3 M© planets in the HZ for more centrally- 
condensed disks, i.e., for larger values of a (not shown 
in Fig. [J). Given the large amount of variation in fZ 
and other uncertainties, we do not consider these limits 
to be firm. However, given the large uncertainties and 
expected variation in / and other parameters, we do not 
consider these limits to be meaningful except in a sta- 
tistical sense. For an ensemble of disks, the fraction of 

> 0.3 M© planets that form decreases significantly for 
low-mass stars. A probabilistic version of the mass limit 
estimate is discussed further in §4. 

3.2. Formation Timescales and Planetary Water 
Contents 

Figure [5] shows the mean formation timescales for HZ 
planets in our simulations. Around 0.2 M© stars, ter- 



restrial planets in the HZ form in a few Myr. This in- 
creases to 20-50 Myr for Sun-like stars, consistent with 
estimates from Hf/W isotopic measurements (e.g., Ja- 
cobsen 2005). The reason for the speedup in accretion 
times in the HZs of low-mass stars is due to a combi- 
nation of the faster orbital speeds in the HZs of low- 
mass stars and the higher surface densities (albeit much 
lower total HZ masses). For instance, if we assume 
L+ oc Mf (a rough fit to Hillenbrand & White 2004), 
then the location of the HZ, rjjz, scales with the stel- 
lar mass as rjjz oc L^J 2 oc M%. The orbital speed in 

the HZ, vhz, scales as vhz oc (M^/thz) 1 ^ 2 oc A/* 1 . 
The surface density in the HZ, "Ehz, scales as S^z oc 
r HZ M * M~ 2a+h oc M" 1 for our case of a = h = 1. 
So, using a very rough approximation that the growth 
time to scales inversely with the product of the orbital 
frequency (~ vhz/thz) and the local surface density 
(Safronov 1969) yields to oc M 7 J 2 for our case. This 
scaling is a very poor match to our simulations, yield- 
ing accretion timescales that are far too short for stellar 
masses below 0.6 M© (dotted line in Fig. [5]). This is be- 
cause the final stage of planetary growth is dominated by 
isolated scattering events between embryos rather than 
accretion from a continuous medium of planetesimals. 
In fact, we notice that a slightly different scaling, with 

to oc {^hzvhz) -1 oc M^ 2 , provides a good empirical fit 
to our simulated formation times (dashed line in Fig. [5]). 
We suspect that this is simply because of the dynam- 
ics of this accretion regime, in which the relevant bodies 
are not on initially crossing orbits, and secular pertur- 
bations are required to excite eccentricities and thus to 
cause orbits to cross. 

Lissauer (2007) argued that the formation time of HZ 
planets around low-mass stars is very short, a few hun- 
dred thousand years for a 1/3 M© star. Our simulations 
confirm the trend that HZ planets form faster around 
low-mass stars. However, Fig. [5] shows that these for- 
mation times do not scale nearly as strongly with stellar 
mass as indicated by Lissauer (he calculated to ~ M®' 2 ). 
The reason for this discrepancy is that Lissauer required 
an Earth-mass planet to form in the HZ, so to accomplish 
this he increased the disk mass by a large factor. Our 
models suggest that a factor of 20-30 increase is needed 
to form a 1 M© planet in the HZ of a 1/3 M© star; as 
expected, that factor corresponds to the approximate dif- 
ference in accretion times between our simulations and 
Lissauer's. However, the fraction of 1/3 M© stars with 
disks massive enough to form a 1 M© planet in the HZ is 
less than 5%, although that fraction increases to about 
15% if a = 3/2 (see Fig. |6] below). In addition, recall 

that the much shallower empirical scaling (to oc M^ 2 
for our case) provides a far better fit to our simulations 
than the simple derivation using the orbital frequency 

(t G oc M; /M for our case; based on Safronov 1969 and 
calculated in similar fashion to Lissauer 2007). 

The mean water content of our HZ planets decreases 
drastically around low-mass stars. In fact, only one out 
of 19 HZ planets formed from a star with M* < 0.6M© 
contained a significant amount of water, compared with 
two out of six HZ planets for M* = 0.8 M© and two 
out of four for M* = 1 M© (see Fig. [2]). The reason 
for this trend - drier planets in the HZs of lower-mass 
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stars - is inefficient dynamical stirring in the low-mass 
disks found around low-mass stars, because more mas- 
sive embryos are needed to increase eccentricities enough 
for significant radial mixing and therefore water delivery 
(Morbidelli et al.2000; Raymond et al. 2004, 2007). Lis- 
sauer (2007) argued that planets in the HZs of M stars 
would be deficient in volatiles for three reasons: 1) col- 
lision velocities are higher; 2) formation timescales are 
faster; and 3) pre-main sequence evolution of the lowest- 
mass M stars is slow, such that the snow line and HZ 
move inward on a Myr to Gyr timescale (see Kennedy 
et al. 2006). We agree with point 1: collision velocities 
are proportional to the orbital speed, so the velocity of 
impactors at infinity (neglecting the escape speed) scales 

with stellar mass roughly as M* ' (see above). We also 
agree qualitatively with point 2: formation timescales are 
indeed likely to be shorter around low-mass stars by per- 
haps an order of magnitude compared with solar-mass 
stars (Fig. [5]). Point 3 is more uncertain: the inward 
movement of the snow line could affect the water con- 
tents of HZ planets if the formation timescale is shorter 
than the snow line's drift timescale. In such a case, ma- 
terial in a given region may be too hot to contain water 
at early times, during accretion. After several Myr, that 
zone may drop below the threshold for water to con- 
dense. If, however, the planets are fully formed by this 
time, then water delivery is impossible. Note that this 
pre-main sequence scenario applies only to very low-mass 
stars. Our results suggest that this is a moot point, as HZ 
planets around such very low-mass stars are unlikely to 
have wide enough feeding zones to accrete water-bearing 
material, even at late times. 

Thus, we argue that terrestrial planets in the HZs of 
low-mass stars are likely to be dry, but for a different rea- 
son than Lissauer. Simply put, very little water-rich ma- 
terial is likely to impact such planets at all, because ra- 
dial mixing is inefficient in low-mass disks preferentially 
found around low-mass stars. The influence of additional 
gas or ice giant planets in the system is not expected to 
increase the water contents of HZ planets, at least not 
from the asteroidal source of water considered here (S. 
Raymond 2007, in preparation). However, migrating gi- 
ant planets (not modeled here) can stir up eccentricities 
and induce the formation of very water-rich planets in 
their wake (Raymond et al. 2006b; Fogg & Nelson 2007; 
Mandcll et al. 2007) . In addition, subsequent delivery of 
water from a cometary source is possible. 

3.3. Effects of Giant Planets 

An obvious criticism of this work is the lack of giant 
planets. Simulations have shown that giant planets play 
an important role in terrestrial planet formation (e.g., 
Wetherill 1996; Levison & Agnor 2003). Giant planets 
on orbits exterior to the terrestrial region, such as Jupiter 
and Saturn, stir up the eccentricities of embryos from the 
outside in via secular and resonant perturbations, while 
mutual scattering of embryos excites eccentricities from 
the inside out (e.g., Raymond et al. 2005b). Accretion 
is suppressed in the vicinity of the giant planet, because 
perturbations between embryos during their growth can 
scatter them into unstable regions such as strong giant 
planet resonances (Wetherill 1994). Indeed, the com- 
bined effects of embryo scattering and perturbations from 



Jupiter and Saturn are thought to be the cause of the 
depletion of the Solar System's asteroid belt (Wetherill 
1992; Chambers & Wetherill 2001; Petit et al. 2001). 

The net effect of giant planets is to increase the ec- 
centricities of embryos during accretion. This can cause 
the width of the feeding zone of terrestrial planets to be 
increased, causing planets to be somewhat more massive 
and less numerous than in the absence of giant planets 
(Levison & Agnor 2003). This increase in planet mass 
can be as high as ~ 30% (or planet masses can be de- 
creased) , but occurs primarily for relatively weak pertur- 
bations, i.e., for less massive giant planets on relatively 
circular orbits (S. Raymond 2007, in preparation). Even 
on circular orbits, external giant planets rarely enhance 
water delivery. Rather, their perturbations clear out ex- 
ternal, water-rich material leaving behind relatively dry 
terrestrial planets (S. Raymond 2007, in preparation). 
Indeed, terrestrial planets in systems with giant plan- 
ets on eccentric orbits are likely to be dry (Chambers 
& Cassen 2002; Raymond et at 2004; Raymond 2006). 
Note that systems with close-in giant planets may also 
contain Earth-like planets, which should be very water- 
rich (Raymond et al. 2006b; Mandell et al. 2007). 

To test the effects of giant planets, we performed eight 
additional simulations including one giant planet on an 
exterior orbit. We chose the 0.6 M Q case for these simu- 
lations, to see if giant planet perturbations could increase 
the mean terrestrial planet mass above 0.3 M©. We used 
the same four sets of embryos described in Table 1 and 
included a giant planet on a circular orbit at 1.3 AU, cor- 
responding roughly to the orbital distance with the same 
temperature (and therefore at a comparable dynamical 
separation from the HZ) as Jupiter for a 0.6 M Q star. 
We ran two sets of simulations, one each for a Neptune- 
and a Jupiter-mass giant planet. 

Planets that formed in the HZ of our 0.6 M Q simu- 
lations without giant planets had masses between 0.06 
and 0.18 M®, with a mean of 0.10 M®. Systems with 
a Jupiter- [or Neptune-] mass giant planet formed, re- 
spectively, HZ planets with masses between 0.05 [0.06] 
and 0.22'[0.10] M®, with a mean of 0.11 [0.09] M®, very 
close to the no giant planets cases. None of the seven 
HZ planets in systems with no giant planets contained 
any water-rich material. Two out of eleven HZ planets 
that formed in the eight giant planet systems contained 
a substantial amount of water. While this is only a small 
fraction of outcomes, it does show that giant planet stir- 
ring can, in some cases, help in water delivery. However, 
in the vast majority of situations, giant planets either 
hinder or have no effect on water delivery (S. Raymond 
2007, in preparation). The formation timescales of HZ 
planets in giant planet simulations were comparable to, 
or even slightly longer than, those with no giant planets. 
In almost all cases, planets with or without giant plan- 
ets reached 75% of their final masses within 10-20 Myr 
(see Fig. [5]). We suspect that the reason giant planets 
do not accelerate accretion in this case is because more 
radial mixing is happening, such that accretion is some- 
what less confined to a given annular zone. This effect 
appears to be small, as less than 20% of planets in the 
giant planet simulations contain water. 

4. DISCUSSION AND CONCLUSIONS 



() 



Raymond, Scalo & Meadows 



Our analysis is based on three key assumptions: 1) ter- 
restrial planets below a given mass are unlikely to sustain 
life on Gyr timescales - following Williams et al. (1998), 
we derive a limiting planet mass of roughly 0.3 M®; 2) a 
protoplanetary disk can be well described by just a few 
parameters (see Eq. 1): the radial surface density slope 
a, the relative disk mass in solids fZ, and the exponent 
characterizing the disk mass-stellar mass relationship, h; 
and 3) the typical planet mass in the habitable zone is 
proportional to the disk mass in that zone (see Eq. 2 and 
Fig. [3]). We performed a suite of twenty dynamical simu- 
lations of the late-stage growth of terrestrial planets that 
confirmed that this scaling does hold for a wide range of 
stellar masses. We assumed in our calculations that no 
giant planets were present. We tested that assumption 
for the 0.6 M© case and found that the effects of exte- 
rior giant planets do not increase the masses or decrease 
the formation times of terrestrial planets, at least for the 
cases considered here. 

Starting from these assumptions, we explored the 
range of parameters that allow planets more massive 
than 0.3 M® to form in the habitable zones (HZs) of 
their host stars. We found that any realistic combina- 
tion of parameters led to a stellar mass - planet mass 
correlation such that HZ planet masses decrease with 
decreasing stellar mass (see Fig. [3]). For a given set of 
parameters (a, h,fZ), one can derive a stellar mass be- 
low which the probability of a habitable planet forming 
in a disk with those parameters decreases to zero (see 
Fig. [4|). For instance, for our fiducial parameters (h = 1, 
a = 1,/ = 1.2, chosen so that the simulations produce a 
mean planet mass of 1 M® for a 1 M® star), this critical 
stellar mass is 0.7 M Q . This would mean that no M stars 
and few K stars should have sustainably habitable plan- 
ets. However, there are some caveats to this approach. 
First, for main sequence stars, it is impossible to go back 
in time to learn the properties of their protoplanetary 
disks, most importantly / and a. Second, for a given set 
of parameters, there still exists a range of potential out- 
comes because of the stochastic nature of the final stages 
of planetary growth (e.g., Wetherill 1996, Raymond et 
al. 2004; see Fig. [2]). Third and most importantly, there 
is no universal set of parameters that applies for all disks. 
A given set of parameter values, corresponding to a spe- 
cific disk or set of disks, has a corresponding planet mass 
distribution. For instance, the above-mentioned 0.7 M® 
limit decreases to 0.31 M® for a more massive disk with 
fZ = 10, and would decrease farther still for ft, < 1 or 
a > 1. Thus, our approach is only viable in a statistical 
sense, by combining likely outcomes from a distribution 
of disks with varying properties. 

What are the values of the crucial parameters? We 
expect h and a to be universal values characterizing the 
distributions of all disks, but / and Z to vary from disk 
to disk. As discussed above, there is considerable uncer- 
tainty in the value of h, which controls the disk mass- 
stellar mass scaling. If h were larger than 2 or smaller 
than 0.5, we expect that the trend would be apparent 
in Fig. 3 in Scholz et al. (2006). Thus, a roughly linear 
disk mass-stellar mass relationship (h = 1) seems likely. 
A range of roughly 2 orders of magnitude in disk mass is 
observed around a given stellar mass (e.g., Eisner & Car- 
penter 2003; Williams et al. 2005; Scholz et al. 2006), im- 
plying a corresponding range in / values. Note, however, 



that there exists a large uncertainty in tying a given / 
value to an absolute disk mass or surface density, because 
of both observational and theoretical uncertainties (e.g., 
Carpenter et al. 2005). There exists a range of 3-5 in 
metallicity Z for thin-disk stars (Nordstrom et al. 2004), 
although the range for currently-forming stars may be 
much smaller, only ^20% (see Cartledge et al. 2006). 
Nearly all disk models and observations suggest that the 
surface density exponent a — 0.5 — 1.5 (e.g., Dullemond 
et al. 2007; Andrews & Williams 2006). Our choice of 
a = 1 as a fiducial case was a compromise between these 
varying estimates. 

Given the large amount of intrinsic variation in the disk 
mass and accretion outcome, as well as the uncertainty 
in h and a, we do not consider any stellar mass limits 
from Fig. [4] to be firm. Rather, our primary result is a 
statistical observation that, in our framework, the frac- 
tion of disks capable of forming planets of a given mass 
in the habitable zone decreases around low-mass stars. 
We investigate a model with fZ distributed as a Gaus- 
sian in log space, with a mean of fZ = 1 and a standard 
deviation of 0.5 dex. For simplicity, all other parameters 
are held fixed, with a = 1, and h—1. Figure [6] shows 
that the fraction of systems that can form a > 0.3 M® 
HZ planet decreases from 80% around solar-type stars to 
42% around 0.7 M® stars to 6% around 0.4 M Q stars to 
only 0.3% around 0.2 M® stars. Disks with more cen- 
trally concentrated mass distributions (larger a values) 
have less difficulty forming massive planets in HZs very 
close to low-mass stars (e.g., Raymond et al. 2005b). For 
our chosen fiducial parameters, we estimate that only 
25% of systems with M* = 0.6 M Q can form planets 
> 0.3 M®, and only about 5% for M* = 0.4 M Q . Disks 
with radial surface density profiles as flat as r -1 / 2 have 
near zero probability for stars with M* < 0.5 M®, while 
probabilities are 2-4 times larger for disks with r~ 3 / 2 
profiles. Reducing the required planetary mass for hab- 
itability to 0.1 Mq increases the probability by a factor 
of 2-4 (see below). 

Our estimate of the minimum planet mass that can 
sustain enough tectonic activity for a period of time that 
we took as 5 Gyr is also very uncertain. It is simple to 
show that if the mass luminosity relation is of the form 
oc Mf, then the derived critical stellar mass M* jCr be- 
low which there should be few planets with masses above 
a critical mass M p i jCr is given by M+ iCr ~ {M p i tCr / f) B , 
with B — 2/[p(2 — a) + 2h]. Thus, a reduction in the 
adopted critical mass for sustained habitability from 0.3 
to 0.2 M® reduces the required value of disk mass nor- 
malization by the same factor, in order to obtain the 
same critical stellar mass. Alternatively, for a given cho- 
sen disk mass normalization, the derived M* lCr varies 
with the assumed planet mass to a power that is between 
about 0.3 to 0.7 for most of the range of parameters de- 
scribed above. The inset of Fig. [5] shows the effect of 
changing our 0.3 M® lower limit, retaining fixed values 
of h = 1 and a = 1. 

A system of three planets with minimum masses be- 
tween 5 and 15 M® was recently discovered orbiting the 
0.31 M star Gliese 581 (Bonfils et al. 2005; Udry et 
al. 2007) . Figures |4] and [6] suggest that such a low-mass 
star is unlikely to form massive terrestrial planets. How 
can we reconcile this? We see two possible explanations: 
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1) the planets of Gl 581 formed in situ from a very mas- 
sive disk, with fZ ~ 30 — 50; or 2) our assumption of in 
situ formation failed in this case: the planets of Gl 581 
formed farther from the star (probably in the icy regions 
of a relatively massive disk, with fZ of at least 5-10) 
and migrated inward (Goldreich & Tremaine 1980). If 
any of the Gl 581 planets were to transit the parent star, 
then it would be possible to derive a rough composition 
(rocky vs. icy: Valencia et al. 2007; Fortney et al. 2007; 
Sotin et al. 2007) , and to distinguish between these two 
formation scenarios (Gaidos et al. 2007; S. Raymond et 
al. 2007, in preparation). Note that, even if the Gl 581 
planets formed in situ, we expect that they would contain 
a substantial amount of water. Following the arguments 
made above, massive disks promote eccentricity growth 
and radial mixing such that planets in the HZ should 
accrete a large amount of water-rich material. 

Even if HZ planets of sufficient mass can form around 
M stars, they face other challenges, including loss of at- 
mosphere by intense stellar activity coupled with their 
small distance to their parent stars, depending on the 
strength of the planetary magnetic field (Lammer et 
al. 2007) , a snow line that takes perhaps a Gyr to get 
close enough to the HZ to give sufficient water (Kennedy 
et aL2006), or loss of volatiles due to larger impact 
speeds and faster formation times (Lissauer 2007; but see 



§3.2 and Fig. [5]). However these latter two points should 
not affect decisions concerning M star planet searches, 
since they only apply to the very lowest-mass stars, 
~ 0.1 — 0.2 Mq, whose apparent brightness and HZ an- 
gular separation are far too small for any planned search 
(Scalo et al. 2007). Unless the masses of disks are larger 
than is currently thought by a significant factor and/or 
the critical planet mass for habitability is considerably 
smaller, only a small fraction of accessible M star sys- 
tems (with Af* > 0.3 - 0.4 M ) should have habitable 
planets that remain habitable for billions of years. 
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TABLE 1 

Initial Conditions for Simulations 1 



Range (AU) Mass(M ffi ) N Timestep (d) HZ (AU) Water line (AU) 



1.0 


0.5-4 


4.95 


75 


6 


0.8-1.5 


2.5 


0.8 


0.25-2 


1.98 


75 


2.5 


0.39-0.74 


1.23 


0.6 


0.12-1 


0.75 


100 


1.0 


0.20-0.37 


0.61 


0,1 


0.06-1 


0.53 


130 


0,1 


0.10-0.19 


0.32 


0.2 


0.03-0.5 


0.13 


190 


0.2 


0.05-0.1 


0.16 



1 Columns are labeled as follows: 'M*' is the stellar mass in Solar masses; 'Range' is the 
initial radial distribution of embryos; 'Mass' represents the total mass in embryos; 'N' is 
the number of embryos; 'Timestep' is the timestep used for the integration in days; 'HZ' 
is the extent of the habitable zone for each case; and 'Water line' indicates the boundary 
beyond which embryos are assumed to contain 5% water by mass. 
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Fig. 1. — Evolution of a simulation around a 0.6 Mq star. Particle size scales with the planetary size as the planet mass 1 / 3 , but is not to 
scale on the x axis. Colors correspond to water contents, from red (dry) to blue (5% water by mass - see color bar). For scale, the planet 
that formed at 0.29 AU in the HZ is 0.06 M® and the planet at 0.42 AU is 0.21 M®. The HZ is shaded in the final panel of the simulation. 
Note that, although the HZ planet is dry, some water delivery did occur, as two water-rich embryos were accreted by the planet at 0.41 
AU. 
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Fig. 2. — Final outcomes of ten simulations, two for each stellar mass chosen, and with the Solar System for scale. As in Fig[T] color 
represents water content, and particle size scales with the planet size, i.e., the planet mass 1 / 3 . The lines under each planet represent the 
radial excursion of a planet over its orbit, i.e. its orbital eccentricity - these values are averaged over the last 10 Myr of each simulation. 
Solar System water contents are from Lodders & Fegley (1998) and Lecuyer et al. (1998), and eccentricities are 3 Myr averages from Quinn 
et al. (1991). The HZ is shaded for each stellar mass. 
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Fig. 3. — Mass of planets formed in the HZ as a function of stellar mass, for a model with h = 1, a = 1, and fZ = 1.2 (so that the 
mean planet mass for a 1 ~Mq star is 1 Mgj). Error bars represent the range of values for HZ planets. The solid curve represents a model 
in which the HZ planet mass scales linearly with the total annular mass in the HZ. The shaded region represents reasonable estimates of 
the limiting planet mass for habitability (0.1-0.5 Mj); our chosen value of 0.3 M^ is indicated with the dashed line. 
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Fig. 4. — Regions of M+ — fZ space in which habitable planets more massive that 0.3 Mg can form, assuming a = 1 and for three 
different values of h. Planets larger than 0.3 Mj) can form above and to the right of each curve. 




Fig. 5. — Formation times of HZ planets in our simulations. Different symbols correspond to the time for a planet to reach a fraction 
(50%, 75%, 90%) of its final mass. Shaded are estimates for the formation time of the Earth, derived from Hf/W isotopic measurements 
(e.g., Jacobsen 2005). The dotted line corresponds to a simple estimate from Safronov (1969), assuming the formation time scales inversely 
with the product of the orbital frequency and the local surface density. The dashed line represents a different, simple model in which the 
formation time scales inversely as the product of the orbital velocity and the local surface density. Both estimates are referenced to 50 Myr 
for 1 Mq. 
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Fig. 6. — Fraction of stars able to form a > 0.3 Mj HZ planet as a function of stellar mass. We assume an fZ distribution that is 
gaussian with a standard deviation of 0.5 dex, and a fixed value of h = 1. The inset shows the effect of varying the limiting mass for 
habitability from 0.1 to 0.5 M®, for a = h = 1. 



